Bose-Einstein source of intermittency in hadronic interactions 



Tadeusz Wibig 
Experimental Physics Dept., University of Lodz, 
ul. Pomorska 149/153, PL-90-236 Lodz, Poland 
(February 1, 2008) 

Abstract 

The multi-particle Bose-Einstein correlations are the source of "intermit- 
tency" in high energy hadronic collisions. The power-law like increase of 
factorial moments with decreasing bin size was obtained by complete event 
weighing technique with gaussian approximation of space-time particle emit- 
ting source shape. The value of source size parameter was found to be higher 
than the common one fitted with the help of the standard Handbury Brown- 

Twiss procedure. 
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The use of intensity interferometry to determine space-time sizes of the particle emitting 
source is a well-established technique of high energy physics. The standard methods based 
on Handbury Brown-Twiss |I| (HBT) effect is to fit the Fourier transform of the source 
space-time density to the two-particle correlation function. The probability of finding one 
of the two emitted particles with the momentum pi and the second with p 2 is given by 

P{i2} = J \^(x 1 ,x 2 ;puP2)\ 2 p(xi) p(x 2 ) d*x 1 d 4 x 2 , (1) 

where x\ and x 2 are the four-positions of the emission points each of them distributed in 
the "source" according to p. If the particles are bosons a symmetrization in the amplitude 
\l/ evaluation leads to the well-known formula: 

P {12} ~ 1 + |.F 12 | 2 , (2) 

where the T\ 2 can be related to the source distribution by 

Tn = J eto'rt p(x) d A x- % = Vl - Pj . (3) 

There could be also other interpretations of T\ 2 . As for an example the one given in Ref. 
0] derived on a basis of the relativistic string fragmentation picture. 

The particular choice of the source space-time distribution (or, more general, the form 
of T\ 2 ) leaves some degrees of freedom here, but the results do not depend very much of 
that choice. The most popular is the gaussian in space and exponential in time emission 
source shape. However for the present work we choose the form of T\ 2 which is known as a 
gaussian parametrisation for its simplicity and because it is Lorentz-invariant: 

^ = C -««W>; Ql = -{p t - Po )\ (4) 

which leads to the well-known formula for the two-particle correlation function 

C 2 (Q 2 ) = 1 +e~^ 2 . (5) 

The Rq can be still interpreted as a measure of space-time extension of the emitting 
source. 
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The idea presented shortly above has extensively been used to analyze different high 
energy physics data since first work |J by Goldhaber et al. Since that time many experi- 
mental and theoretical efforts have been made. The different source shapes were examined, 
some fine effects were predicted. Some difficulties were also found in the interpretation of 
the source shape while the source is moving very fast with respect to the laboratory system. 
However the main idea of HBT effect remains unchanged. 

In the mid-eighties due to the work of Bialas and Peschanski Q a new interest for the 
particle correlation has arisen. The phenomenon called "intermittency" was found in the 
very small phase-space bin size analysis. Since the first measurements the experimentally 
available smallest bin size is reduced more than order of magnitude but, what is even more 
important, new techniques to study fine structures were developed. The "intermittency" of 
the particle creation process seen by Bialas and Peschanski, which is in fact the fractal (self- 
scaled) behaviour of the multiparticle correlation measures at the very small phase-space 
scales, contradicts the standard Bose-Einstein statistics driven description given by Eqs. (0) 
and (|3p. The intermittent picture of hadronic creation was also inconsistent with existing 
models of particle production (like e.g. LUND hadronization model). The intermittent 
models like a-model j|,|5||, geometrical branching model ]6[], one- dimensional model of inter- 
mittency by Dias de Deus were invented but none of them achieved such a completeness 
and predictivity as high energy physics standards (LUND or DPM-type models). On the 
other hand, the treatment of the "intermittency" as a real new phenomenon was still not 
so obvious. In Ref. || different data sets were examined and as the last conclusion it is 
stated that the intermittency is caused by Bose-Einstein correlations in addition to a mech- 
anism responsible for the power-law behaviour, in Ref. || authors claim that the observed 
"intermittent-like" behaviour of moments of multiplicity distributions can be understood as 
an effect driven by quantum statistical properties of the particle emitting system and it does 



not necessarily imply evidence for intermittency. The title of the Ref. |]T0| "Has intermit- 
tency been observed in multi-particle production?" is, to some extend, a good question still 
nowadays. 



As it has been said, the existing data gives a possibility to study intermittent (power- 
low) behaviour of factorial moments in more than two decades wide phase-space distance 
measure (however it will be defined: rapidity, momentum or four-momentum difference, 
box volume etc.). It is clear now that the classical picture describing the HBT effect (Eq. 
(0)) is not valid. The factorial moment analysis were performed for e + e~ annihilation into 
hadrons experiments, for hadronic collisions 7r + /K + p at y/s = 22 GeV and at much higher 
energy \/s = 630 GeV pp. There are also data from pM and MM experiments. All they 
shows more or less definite a power-law like dependence on the bin size. However there is 
also very clear signal about the like and unlike charge difference of the correlation strength 
which suggests its Bose-Einstein origin. The possibility to achieve an agreement between 
those two, on the first sight contradicting, experimental facts will be discussed in the frame 
of common quantum physics. 

It should be remembered that Eq. @ was obtained in the case when only two particles 
were emitted from the source. That situation is of course different when one has got to 
do with multi-particle source [ XT | . In some particular cases (when there are really a small 



number of particles emitted in the large phase-space volume) the two-particle correlator 
given by Eq. (^) still can be used as at least a first approximation. But when one wants to 
look closely at the high multiplicity events or to study multiparticle correlations Eq. @) has 
to be modified. 

When n identical bosons are emitted the probability of the particular momenta configu- 
ration {p;} is given by: 

~ X^Ml) -^2(7(2) ••• Fna(n), (6) 

where a is a permutation of a sequence {1,2, ...,n}, cr(i) is the i-th element of this permu- 
tation and the sum is over all n! permutations. 

To see what real difference is introduced by such complete treatment the two-particle 
correlation function like that in the Eq. (0) in the case of three particle emitting source is 
written explicit below: 



P 3 ~ 1 + + l-^isl + l^2 3 | + 2 l^ia | |^i 3 | |^2 3 |- (7) 

If all three particles are very close to each other the statistical weight of such events tends 
to 6 = n!. The limit for two-particle correlator in n-particle emitting source is n! not 2 like 
it comes form Eq. (@). The same limit was obtained in Ref. [Q] but it was interpreted as a 
limit of n-th factorial moment. The multiplicity distribution in the very small phase-space 
bin tends to the geometrical one which, on the other side, can be treated as Bose-Einstein 
statistics driven multiplicity distribution while n — > oo, 5 — > with nS = const. 



Quite different approach to Bose-Einstein phenomenon is discussed in Ref. ||12|| . The 
authors argued for the local nature of the Bose-Einstein effect. In general, their treatments 
leads to the weighing procedure with the event probability proportional to: 

Pin}- E ( 1 + inf)- (8) 
all pairs 

The definition of T' in Eq. @ is not given by Eq. (f|) but is based on string fragmentation 
picture. However the difference is rather in the physical interpretation than in the general 
behaviour. It should be noted that Eq. @ overestimates the very close particle limit. 
It is there equal to 2 n(n ~ 1 * ) / 2 . The arguments for such a treatment are discussed in Ref. 
(similar attempt is presented in Ref. |13|| ) and will not be discussed here. One of 



the arguments not given there but of the practical importance is that the above idea can 
be easily incorporated into the Monte-Carlo event generator. It was in fact done in the 
LUBOEI subroutine which is a part of LUND hadronization scheme JETSET 7.3. The 
general difference between the Ref. |l2j strategy and proposed in the present paper is in 
the fact that sum in the Eq. @ is performed over permutations of the particle ensemble 
in which only two particles are exchanged (locality of Bose-Einstein interaction) while in 
our treatment all event permutations can give a contribution to the event weight (global 
Bose-Einstein approach). The importance of many particle exchange contributions will be 
discussed later on. 

The problem with complete weighing procedure is also a practical one. The sum over n! 
elements can be performed easily for about of ten particles or less. For higher multiplicities 



the calculation time rise tremendously. But it is quite clear than for the two very distant 
particle exchange the contribution coming from all permutations concerning that particular 
exchange is negligible. The algorithm was invented to omit all the negligible permutations 
and calculations of the weights according to Eq. (|6]) became possible also for larger multi- 
plicities. In the present paper only the data from NA22 experiment will be analyzed. The 
mean charged particle multiplicity is of order of 8 and the largest like type boson multiplicity 
(in one chain, as will be discussed later) in the sample of about 500000 our Monte-Carlo 
generated events does not exceed 15. 

To study the influence of Bose-Einstein weighing method on the shape of the fluctuations 
in small bins the sample of events in the "world of absence of Bose-Einstein correlation" is 
needed. There is a number of Monte-Carlo generators which can be used to get this. In 
the present work the one called Geometrical Two-Chain was used. It is described in details 



in Refs. [fT3 |. The advantage of that generator is the minimum of correlations introduced 
there. The ones existing are due to the conservation requirements (charge, barion number, 
strangeness, momentum and energy), the resonance production and the large scale clustering 
due to chain mass distribution in the model. There are also correlations connected with 
the hadronization procedure adapted: the transverse momentum is conserved locally in 
the fragmenting chains so the subsequent hadrons incline to have the negatively correlated 
momenta perpendicular to the interaction axis. Our chain fragmentation picture leads also 
to ordering in rapidity of subsequently produced hadrons. All that features are present in 
most of the models working on the partonic level. The last but very well seen specially for 
large bin sizes is a contribution related to non-poissonian multiplicity distribution in the 
multiparticle production. 

The main interaction characteristics are very well reproducible by the generator as it was 



shown in Refs. [ |T4|j . 

About 500000 of non-single-diffractive events for 7r + and K + interactions with proton 
at laboratory momentum of 250 GeV/c were generated and combined to get the reference 
sample without Bose-Einstein correlations included. Then for each event the weight was 
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calculated according to Eqs. (f|) and (|6]). In principle the Bose-Einstein weighing procedure 
could change the multiplicity distribution (what was one of the argument against global 
treatment of Bose-Einstein correlation in Ref. [12| ). To avoid this the weights were renor- 



malized to get the average value of the weights for n identical bosons equal to 1 and these 
were used afterward. The detail comparison with the experimental data leads to the conclu- 
sion that if the Bose-Einstein symmetrization were performed for the whole events then the 
correlations are too strong for very small bin sizes. In our model there is only one parameter 
to be adjusted, correlation radius Rq, while in the standard HBT procedure there is also 
the incoherence parameter which allows to make softer the correlation strength. In the Ge- 
ometrical Two-Chain model particles are produced by the fragmentation of two well-defined 
chains so there is a natural subdivision of all secondary particles to two distinct classes. 
To make the correlation weaker there is a possibility to symmetrize amplitudes \1/ not over 
all particle exchanges but only over the exchanges of the particles produced from the same 
chain. 

The very convenient variable to study the two-particle correlation is the differential form 
of the second factorial moment as it was used in Ref. ||15||. The definition using density 



integral method |T6 



is: 



D 2 (Q 2 ) = j±- 2 £ 6(Q 2 - Q%) x 0(QJ - Q 2 + 5), (9) 

Kj 

were B is the Heaviside unit step-function and Norm is a normalization term defined by the 
so-called "mixed events" technique. The particles used for the normalization were chosen 
randomly from the all event "pool" of the large number of generated interactions, ensuring 
that they belongs to different real events. To avoid in the reference sample the correlations 
due to non-poissonian multiplicity distribution in hadronic interactions the multiplicity in 
the mixed events were taken from Poisson distribution with the average value the same as 
in M-C generated events. 

In the particular NA22 experimental data which we want to compare with the rapidity 
cut |y| < 2 has been used. Thus in all the calculations the same cut is applied. In the 



experimental procedure there was also not possible, in general, to determined the particle 
masses so all the particles (except low energy proton and very energetic particles in K + 
induced interactions pi a b > 150 GeV/c) were treated as pion. The same procedure has been 
used in our analysis of the Monte-Carlo events. The experimental accuracy of particle four 
momentum difference determination described in Ref. |T7] was taken into account in the 



calculations as well. The calculations od D2 was performed for all charged particles as well 
as for like and unlike charge combinations. The results are presented in Fig. [I] by the solid 
line. The remained correlations produced in the Geometrical Two-Chains model, which were 
indicated above, leads to the outcome depicted by the dotted line. It represents the result 
of the correlation calculations without Bose-Einstein weighing. 

It is seen that the power-low-like behaviour of D2(q 2 ) is quite well reproduced by 
our weighing method. The small overestimation of the unlike particle correlator at four- 
momentum difference of about Q 2 ~ 10 -2 -j- 10 _1 (GeV/c) 2 is a consequence of the strict 
ordering in rapidity of the chain fragmentation products which introduces always between 
close (in rapidity) like type charged hadron the one with the opposite sign. The four- 
momentum difference of that unlike charged pairs is determined by transverse momentum 
distributions so the effect does not influence the very small bin size analysis. 

However, the main argument for intermittency comes from the analysis of the higher 
multiplicity correlation measurements. To study this effects the correlation measures have 
to be defined for three- and more particle systems. The most commonly used variables 
are the factorial moments. For practical purposes the best method of factorial moment 
calculations is again the one proposed in Ref. density integral method. 



W) = T^-e< E n e(g 2 -^ 1)l(fc2) ), (io) 

l\Ullll i(l)<i(2)<...<i(g) all pairs (i(kl) ,i(k2)) 

with the normalization by mixed event technique again. Results of our calculations are 
presented in Fig. ^. The power-law like increase of factorial moments with decreasing bin 
size is again quite well reproduced in the whole range of Q 2 measured experimentally. 

In Figs, m and ^| the results of event weighing defined by Eqs. (f|) and (||) with the sum 
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over permutations with only one particle pair exchange are also presented (by the dashed 
line). As it is seen the effect (for the same value of Rq parameter in Eq. (^) ) is much 
weaker. This illustrate the importance of global treatment of Bose-Einstein correlation. 
The introduction to the sum of the weights very many relatively small terms leads to really 
great increase of the effect. 

To reproduce the shape of D 2 (Q 2 ) and F 2 (Q 2 ) dependencies measured by NA22 exper- 
iment the value of the parameter Rq in Eq. @) had to be adjusted. The large statistical 
fluctuations of the weights influence the estimation of source size parameter so the accuracy 



achieved is not higher that 10%. In Ref. |15j the source size was found using the standard 
technique of HBT effect (Eq. (||) ). The value found there was (0.82 ± 0.02) fm. Our com- 
plete weighing procedure gives stronger correlations ( even after weights renormalization ) 
so the value of R used to obtain the results given in Figs, p] and || is about 50% higher 
what gives the source radius of about 1.25 fm in the gaussian approximation (Eq. (f|) ) 
interpretation. 

To summarize, the importance of the global treatment of Bose-Einstein correlation has 
been shown. The symmetrization over all permutations leads to the power-low like behaviour 
of factorial moments in the four-momentum difference regions where they are measured 
experimentally. The more detail analysis is in progress and the results will be presented 
elsewhere. 
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FIGURES 

FIG. 1. The differential second factorial moments for a) all charged, b) like-charged and c) 
unlike-charged pairs as a function of four-momentum difference. The data points are from NA22 
experiment. Solid line represents the result of our complete Bose-Einstein weight method, dotted 
shows the correlations in used Geometrical Two-Chain model. Dashed line is for the sum in Eq. 
(6) over only one pair of boson exchanges. 

FIG. 2. The factorial moments a) for negatives and b) for all charged particles as a function 
of four-momentum difference. The data points are from NA22 experiment. 
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